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Vibrational exciton-mediated quantum state transfert: a simple model. 
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A communication protocol is proposed in which quantum state transfer is mediated by a vibra- 
tional exciton. We consider two distant molecular groups grafted on the sides of a lattice. These 
groups behave as two quantum computers where the information in encoded and received. The 
lattice plays the role of a communication channel along which the exciton propagates and interacts 
with a phonon bath. Special attention is paid for describing the system involving an exciton dressed 
by a single phonon mode. The Hamiltonian is thus solved exactly so that the relevance of the 
perturbation theory is checked. Within the nonadiabatic weak-coupling limit, it is shown that the 
system supports three quasi-degenerate states that define the relevant paths followed by the exciton 
to tunnel between the computers. When the model parameters are judiciously chosen, construc- 
tive interferences take place between these paths. Phonon-induced decoherence is minimized and a 
high-fidelity quantum state transfer occurs over a broad temperature range. 

PACS numbers: 03.65.Yz,03.67.-a,63.22.-m,71.35.-y 



I. INTRODUCTION 

Over short length scales, high-fidelity quantum-state 
transfer (QST) from one region to another is a funda- 
mental task in quantum information processing 1 . It is 
required to ensure a perfect communication between the 
different parts of a quantum computer (QC) or between 
adjacent QCs. In that case, solid-state based system is 
the ideal candidate for scalable computing for, at least, 
two main reasons 2 -. First, in a lattice, the information is 
encoded on elementary excitations that naturally propa- 
gate owing to inherent interactions between neighboring 
sites. QST is thus spontaneously achieved without any 
external control. Second, no interfacing is needed be- 
tween the QCs and the communication channel (CC) that 
involves the same elementary excitations. Because the 
CC depends on the way the information is implemented, 
different strategies have been elaborated, such as optical 
lattices 3 -, arrays of quantum dots^, conducting polymers 5 - 
and spin networks^—. However, it has been shown that 
qubits may be encoded on intramolecular vibrational 7 - - — 
so that vibration-mediated QST is a promising alterna- 
tive for quantum information processin g 24 ' 25 . 

In this paper, we propose a protocol in which a molec- 
ular lattice plays the role of the CC, QST being me- 
diated by high-frequency vibrational excitons. Indeed, 
molecular lattices exhibit regularly distributed atomic 
subunits. Owing to dipole-dipole interactions, the energy 
of a specific internal vibration delocalizes between these 
subunits giving rise to narrow-band excitons. Although 
many properties of the excitons have been studied 2 ^— , 
their potential interest for QST has been suggested very 
recently-^ 7 - - — . The main idea is to encode the informa- 
tion on a vibrational qubit defined as a superimposition 
involving the vacuum and one-exciton states. 

Unfortunately, the exciton does not propagate freely 
along the lattice. It interacts with a phonon bath that 
tends to destroy the coherent nature of any vibrational 



qubit. The exciton properties are thus described by a re- 
duced density matrix (RDM) whose behavior is governed 
by a generalized master equation (GME) 50 ~ 52 . To under- 
stand quantum decoherence, we have studied excitonic 
coherences, i.e. the RDM matrix elements that measure 
the ability of the exciton to be in a qubit state^r— . The 
dynamics was described using a Frdhlich model 57 within 
the nonadiabatic weak-coupling limit, i.e. a common sit- 
uation for vibrational excitons. It has been shown that 
high-fidelity QST requires the use of a finite-size system. 

Indeed, in an infinite lattic o 53 ' 54 , the phonons behave 
as a reservoir. The Markov limit is reached and stan- 
dard GME approaches can be used. Dephasing limited- 
coherent motion takes place so that the coherences local- 
ize preventing any efficient QST. By contrast, in a finite- 
size lattice, a strong non-markovian regime occurs result- 
ing in the breakdown of GME method s 55 ' 56 . To overcome 
this problem, perturbation theory (PT) has been applied 
successfully^ 7 - - — . It has been shown that the phonons 
evolve differently depending on whether the exciton oc- 
cupies the vacuum or an excited state. Exciton-phonon 
entanglement takes place resulting in the decay of the 
excitonic coherences. Nevertheless, the decoherence de- 
pends on the nature of the one-exciton state. 

In particular, for odd lattice sizes, the coherence in- 
volving the state located at the band center survives over 
an extremely long-time scale, even at high temperature. 
Unfortunately, this state is a stationary wave unable to 
carry information between the lattice sides. Nevertheless, 
we can take advantage of its robustness against decoher- 
ence to define a protocol quite similar to that introduced 
by Yao et al. in spin chains^ 6 -. The main idea is to con- 
sider QST between two distant molecular groups that 
are grafted on the lattice sides. The structure must be 
designed so that a vibration of each molecular group is 
resonant with the robust state of the lattice while remain- 
ing insensitive to the phonons. The grafted groups play 
the role of the QCs whereas the lattice defines the CC. 
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FIG. 1: Communication protocol between two distant com- 
puters x — and x = L. The information is carried by an 
exciton that propagates along a finite-size lattice. 



Dipole-dipole interactions couple the QCs with the CC 
and allow the exciton propagation so that neither control 
nor interfacing is required. 

In such a confined system, strong non-markovian ef- 
fects will occur so that the dynamics will be addressed 
within the effective Hamiltonian concept provided by 
P T 58 i 59 . However, the fundamental question arises 
whether PT is relevant or not, depending on the model 
parameters. To answer that question, we shall restrict 
our attention to a system in which only the lowest fre- 
quency phonon mode (LFPM) is considered, this mode 
being responsible for the largest perturbation of the exci- 
ton dynamic o 47 ' 48 . Due to its simplicity, this model can 
be solved exactly so that the PT performance can been 
checked easily. 

The paper is organized as follows. In Sec. II, the 
exciton-phonon Hamiltonian is described and QST is for- 
mulated in terms of the excitonic coherences. Then, PT 
is applied to evaluate the system eigenproperties and 
to determine an approximate expression of the exciton 
RDM. In Sec. Ill, a numerical analysis is performed in 
which PT is compared with exact calculations. The re- 
sults are discussed and interpreted in Sec. IV. 



II. THEORETICAL BACKGROUND 

A. Model Hamiltonian 

We consider a CC formed by a ID molecular lattice 
that contains an odd number of sites N (Fig. 1). Each 
site x = 1, N is occupied by an atomic subunit whose 
internal dynamics is described by a two-level system. Let 
\x) be the first excited state of the xth two-level system 
and wo the corresponding Bohr frequency. The vacuum 
state |0 CC ) describes all the two-level systems in their 
ground state. In the CC, the exciton Hamiltonian is de- 
fined in terms of the hopping constant <£> as (H = 1) 

N N-l 

H oc = ^2u \x){x\ + ^2mx + l){x\ + \x){x + l\). (1) 

X— 1 X— 1 



Owing to the confinement, one-exciton states are N sta- 
tionary waves with quantized wave vectors = kir/L, 
with k = 1, .., N and L = N + 1, expressed as 

N n 

l^ fc ) = E vi sin(iffex)|x) - (2) 

x=l 

The corresponding eigenenergies lu® = o; + 2$ cos(Kk) 
form N discrete energy levels centered on u>q, i.e. the 
energy of the so-called robust stationary wave \*-Pl/2)- 
The phonons of the CC refer to the external motions 
of the lattice sites that behave as point masses M con- 
nected via force constants W . They define N normal 
modes with wave vectors Q p — pit / L and frequencies 
Q, p = fl c sm(Q p /2), withp= 1,..,N and n c = y/4W/M. 
Restricting our attention to the LFPM p = 1, the phonon 
Hamiltonian is defined as H B = fta^a, with ft = Q,\, (V 
and a being standard boson operators. In the phonon 
Hilbert space £b, the eigenstates are thus the well-known 
number states \n). As detailed previously^ 7 -, the exciton- 
phonon interaction V = M{a* + a) favors the exciton 
scattering from \cpk) to |v?k±i) mediated by phonon ex- 
changes. The matrix elements of the M operator are 
defined as 

Mfcfc/ = 7?(£fc,fc'+i + 4,fc'-i), (3) 

where 77 = [(E B Q./L)(l - {Q/n c ) 2 )] 1/2 , E B being the 
small polaron binding energy expressed in terms of the 
coupling strength \ as Eb = X 2 /W ■ 

The QCs are formed by two molecular groups x = 
and x = L whose internal dynamics is described by a 
two-level system (Fig. 1). Let |0) (resp. \L)) be the first 
excited state of the Oth (resp. Lth) molecular group and 
ujs the corresponding Bohr frequency. The vacuum state 
\(Z) qc ) describes the two groups in their ground state. Lo- 
cated far enough from each other, these groups do not 
interact. Their internal dynamics is governed by the 
Hamiltonian H qc = w s (|0}(0| + \L)(L\). 

The QCs, assumed to be sufficiently far from the lat- 
tice, are insensitive to the phonons. By contrast, owing 
to dipole-dipole interactions, they interact with the lat- 
tice exciton. This interaction originates in a vibrational 
energy transfer between the group (resp. L) and the 
lattice site x = 1 (resp. x = N) as 

W = <M|0)(1| + |1><0] + \L)(N\ + \N)(L\), (4) 

where $5 = e$, e being the strength of the coupling 
between the QCs and the CC. 

Within this model, the exciton dynamics is governed 
by the Hamiltonian Ha = H cc + H qc + W. It acts in the 
Hilbert space £a whose dimension is N + 3. £a is gen- 
erated by the vacuum |0) = |0 CC ) |0gc) an d by N + 2 
one-exciton eigenstates associated to the eigenfre- 

quencies w M , with fi — 0, L. The exciton Hamiltonian 

is thus rewritten as Ha — X^=o a -vl' ! /v)(' i /vl- To deter- 
mine these one-exciton states, we first fix to luq so 
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FIG. 2: Exciton energy spectrum for e — (left side) and for 
e / (right side). 



that a resonance occurs between the two localized states 
|0) and \L) and the robust stationary wave \^Pl/2) (Fig- 
2). Then, e is assumed to be sufficiently small so that 
off-resonant interactions between the localized states and 
the remaining stationary waves can be disregarded. This 
condition is satisfied provided that e <C ity/2/L. As a 
result, the one-exciton eigenstates are split into two 
groups (Fig. 2). 

The first group contains N — I states that reduce 
to the stationary waves \<fk), with k ^ L/2. For 
fi = 0, ...,L/2 — 2, they are defined as |Vv) = IV/x+i) 



whereas for pL = L/2 + 2, L, they reduce 



to \ipfj_) — l^/i— l) = w u-i)- The second group results 
from the hybridization between the degenerate states |0), 
\ ( Pl/2) an d 1-^)- Indeed, |0) (resp. \L)) interacts with 
\ ( Pl/2) through the coupling constant g = e$>y/2/L (resp. 
g' = g sm(Nir/2)). When e =^ 0, these couplings raise the 
degeneracy so that the eigenstates arc defined as 
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|0) 



V2 V2 



\L), 



(5) 



where Ajv = sin(_/V7r/2) and with the convention /1 = 
L/2 ^ 1 = ± and fi = L/2 = o. The corresponding 
eigenfrequencies are expressed as 



(6) 



ui± = uo±2e$/VL 

iO a = UJ . 



\ip±) and \ip ) define three quasi-degenerate states that 
involve a superimposition of the localized states, either 
symmetric or antisymmetric. Note that \ip ) is exactly 
located at the band center and it does not depend on 
the stationary wave \ipL/2)- This is no longer the case 
for |t/> + ) and \ip—) that lie just above and just below the 
band center, respectively. 



The exciton-phonon Hamiltonian is thus written as 
H = H + V where H = H \ + Hg is the unperturbed 
Hamiltonian. It acts in the Hilbert space £ = £a <8> £b 
that is partitioned into independent subspaces as £ — 
£0 © £\. In the zero-exciton subspace £0, V — so 
that the unperturbed states are eigenstates of H. They 
correspond to tensor products |0) <g> \n) that describe n 
free phonons. In the one-exciton subspace £\, the unper- 
turbed states l*"^) = IVvJ ® \ n ) refer to free phonons 
accompanied by an exciton in state Because V 

turns on, they are no longer eigenstates of H. The ex- 
act eigenstates with eigenenergies E i: correspond to 
entangled exciton-phonon states. 

The coupling V favors exciton scattering from j^) to 
through phonon exchanges. The allowed transi- 
tions are specified by the selection rules (ip^M |Vy) 7^ 
(Eq.Q). Because [ip^) interacts with the phonons 
through its dependence with respect to the stationary 
waves, these rules are quite similar to those that charac- 
terize the CC (|Vv) i s coupled with |Vv±i))- However, 
two main differences occur. First, because \ip ) only de- 
pends on the localized states, it remains insensitive to 
the phonon bath. Second, a splitting occurs for the tran- 
sitions involving \i/j±). Indeed, because \4>l/2±2) reduces 
to \<Pl/2±x)i it interacts with and | ?/>-), both states 
depending on \<Pl/2)- Nevertheless, within the nonadia- 
batic limit, i.e. for 4$ < fl c and e <C ny/2/L, the allowed 
transitions do not conserve the energy. There is no res- 
onance between the coupled unperturbed states so that 
second order PT can be applied to treat V in the weak 
coupling limit (Eg <C 3>). As detailed previousl y 47 ' 48 , 
PT is valid at temperature T provided that L < L* with 
L* m 0.20,1/ E B k B T. When L > L*, quasi-resonances 
take place between the unperturbed coupled states and 
PT breaks down. 



B. Excitonic coherences and QST 

Without any perturbation, the CC and the QCs are in 
thermal equilibrium at temperature T. Assuming that 
loq ^> fcflT, all the two-level systems are in their ground 
state. By contrast, the phonons form a thermal bath de- 
scribed by the density matrix ps = exp(— /3Hb)/Zb, Zb 
being the phonon partition function (/3 = l/k B T). In 
that case, one assumes that the internal vibrations inter- 
act with an external source so that the exciton is initially 
prepared in a state \iPa) 7^ |0)- This step is supposed to 
be rather fast when compared with the typical time evo- 
lution of the phonons. The full system is thus brought in 
a configuration out of equilibrium and its initial density 
matrix is defined as p(0) = \iPa)(" > Pa\ ® Pb- 

To study QST between the QCs, \^>a) defines a qubit 
implemented on the molecular group x = as 



h/u)=a|0)+/?|O), 



(7) 



where \a\ 2 + |/3| 2 = f. Our aim is thus to measure the 
ability of the system to freely evolve in time so that this 
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initial qubit is copied on the second QC x — L (Fig. 1). 
Whatever its duration, QST must be realized with the 
largest fidelity despite the coupling with the phonons. To 
define this fidelity measure, different objects have been 
introduced^, one of the most widely used being certainly 
the so-called average Schumacher's fidelity^. Here, we 
restrict our attention to the excitonic coherences. In- 
deed, the exciton properties are encoded in the RDM 
a{t) = Trs[p(t)], where Ttb is a partial trace over the 
phonon degrees of freedom. The coherences are thus the 
off-diagonal matrix elements a X 0(t). They provide infor- 
mation about the ability of the xth two-level system to 
develop a superimposition between its ground state and 
its first excited state a time t. 

When the initial qubit is implemented on the Oth QC, 
the coherence 000(0) is turned on at time t = 0. There- 
fore, the ability of the Lth QC to develop a superim- 
position involving |0) and \L) at time t is given by 
0X0 (*) = G LO {t)a O0 (O) with 



G L0 (t) = (L\Tr B [ PB e lHBt e- lHt ] |0). 



(8) 



The effective exciton propagator GLo{t) generalizes the 
concept of transition amplitude. It yields the probability 
amplitude to observe the exciton in \L) at time t given 
that it was in |0) at t — 0. Its effective nature results 
from the fact that the exciton interacts with the phonons 
during its transition. 

The effective propagator is the central object of the 
present study. The condition |G_lo(OI — 1 reveals that 
the Lth QC reaches a state at time t that is equivalent to 
the initial state, to a phase factor. Note that this condi- 
tion is exactly the QST fidelity when the coupling with 
the phonons is disregarded?.. Consequently, depending 
on the value of the model parameters, studying the max- 
imum value of |Gz,o(^)| provides key information about 
the fidelity of the QST. 



C. Perturbation theory 

In its operatorial formulation^, standard PT involves 
a unitary transformation that provides a new point of 
view in which the exciton-phonon dynamics is described 
by an effective Hamiltonian. The key point is that this 
Hamiltonian is diagonal in the unperturbed basis. Quite 
powerful to treat finite-size lattice a 47 ' 48 , this approach 
breaks down in the present situation because the unper- 
turbed Hamiltonian exhibits quasi-degenerate states. For 
small e values, n ) and \^°_ „) have almost the same 
energy. Although they do not directly interact through 
V, they are coupled with the same unperturbed states. 
Consequently, effective couplings occur between quasi- 
degenerate states resulting in errors in the calculations 
of the corrected energies. 

To overcome this problem, quasi-degenerate PT is 
applied^ (Appendix A). To proceed, we take advan- 
tage of the fact that the effective couplings conserve the 
phonon number. Therefore, our procedure involves a 



transformation U = exp(S') that generates a new point 
of view in which the effective Hamiltonian H = U HU' is 
block-diagonal in the unperturbed basis. The generator 
S is expanded as a Taylor series in the coupling V so that 
H becomes diagonal with respect to the phonon number 
states, only. Up to second order, it is written as 



H = H A + SH A + + Stya^a 



(9) 



where SHa and <5f2 are operators in Ea whose matrix 
elements are defined as (in the unperturbed basis 



1^ M^M^ 



SO. 



1 x M^^M^! 



M^M^ 2 



+ 



1 \ ~< Mj^Mjjjj^ 



2 ^— ^ uj n , — LU n 



9}' 



SHa is the correction of the exciton Hamiltonian ow- 
ing to the coupling with the phonons. It results from 
the spontaneous emission of a phonon during which the 
exciton realizes a transition from {ip^} to |Vv)- How- 
ever, in the nonadiabatic limit, the energy is not con- 
served during the transition. The emitted phonon is im- 
mediately reabsorbed and the exciton realizes a second 
transition from \ip^} to \ip^ 2 }- In other words the exci- 
ton does no longer propagate freely but it is dressed by 
a virtual phonon cloud. This dressing renormalizes the 
excitonic energies w M by an amount Suj^ — SHa^^- In ad- 
dition, it induces effective interactions SHa^^ between 
distinct excitonic states that can no longer be neglected 
for quasi-degenerate states. 

Similarly, SQ defines the correction of the phonon fre- 
quency. It has two origins. First, the phonon can be 
absorbed giving rise to excitonic transitions. Because 
this process does not conserve the energy, the phonon is 
immediately re-emitted. Second, the phonon can induce 
the stimulated emission of a second phonon during which 
the exciton realizes transitions. But, as previously, the 
emitted phonon is immediately reabsorbed. Both mech- 
anisms are virtual processes indicating that the phonons 
are dressed by virtual excitonic transitions. 

To diagonalize H for each phonon number, we use the 
fact that 5£l is smaller than SHa within the weak cou- 
pling limit. Consequently, let |x„) be the eigenstates of 
Ha + SHa and Cj v the corresponding eigenfrequencies 
(v = 0,...,L). Up to second order in V, one introduces 
<5il„ = (xv\5Q\Xv) as the correction of the phonon fre- 
quency induced by the exciton that occupies the state 
\Xu)- Within these notations, H is rewritten as 



H 



L 

E 

v=0 



&u\Xv){Xv \ + H 



(11) 
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where = (£1 + Sfl^a^a is the Hamiltonian that gov- 
erns the phonon dynamics when the exciton is in the state 
\Xu). 

In the new point of view, the exciton-phonon dynam- 
ics is thus governed by the effective Hamiltonian H that 
is diagonal in the basis \xv) ®\n). Its eigenvalues de- 
fine the system eigenfrequencies up to second order in 
V as E^n — LJ y + n(Q + 5Sl v ), H does no longer refer 
to independent excitations but it characterizes entangled 
exciton-phonon states. A state \%v) describes an exciton 
dressed by a virtual phonon cloud whereas the number 
state \n) describes n phonons clothed by virtual excitonic 
transitions. This entanglement is clearly highlighted in 
the starting point of view in which the eigenstates do no 
longer factorize as \^ v , n ) — U^\xv) ® 

As detailed previously^, PT is particularly suitable for 
deriving an approximate expression for Gxo(i)- To pro- 
ceed, we first introduce U and diagonalize H in Eq.©. 
Then, we define the effective density matrix 



exp [-(f3n + iSn v t)a<a\, (12) 



where Z^\t) = (1 - exp [-(/3fi + i8Q v t)])- x . Strictly 

speaking, p^ (t) is not a density matrix since it yields 
complex values for the phonon population. However, 
it is isomorphic to ps with the correspondence /3f2 — >• 
/3f2 + i5Q u t and it provides averages equivalent to ther- 
mal averages. Consequently, after simple algebraic ma- 
nipulations, Gio(i) is rewritten as^ 



Gw(t) = J2 



Z B 



exp[— iCj v t] x 



(13) 



Tr } 



[ P %\t){L\Ut(t)\Xu)(Xu\U v (0)\Q) 
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FIG. 3: Ei - E v , n vs Ei for e = 0.01, x = 10 pN and L = 10. 
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where U v (t) — e b l JJe~ l b *. Expanding U in a Tay- 
lor series with respect to V, one finally obtains the second 
order expression of G^oWj as detailed in Appendix B. 



III. NUMERICAL RESULTS 

In this section, numerical calculations are carried out 
to show the relevance of PT for describing QST. To pro- 
ceed, the previous formalism is applied to amide-I vi- 
brations in a-helices, a system for which the parameters 
are well-known2&2& : uj = 1660 cm" 1 , W = 15 Nm" 1 , 
M = 1.8 x 1CT 25 kg, fl c = 96.86 cm' 1 and $ = 7.8 cm" 1 . 
To avoid PT breakdown, the size is set to L = 10 and 
the exciton-phonon coupling strength is fixed to % = 10 
pN. Special attention will be paid for characterizing the 
influence of the coupling e between the QCs and the CC. 

The difference AE = Ei — E u _ n between exact and ap- 
proximate energies is shown in Fig. 3. For e = 0.01, PT 
provides a very good estimate of the energy spectrum 
over a broad energy scale. The smaller the energy is, the 



better is the agreement. Of course, AE increases with 
Ei indicating that the PT accuracy decreases with the 
phonon number n because V scales as i/n. For instance, 
for Ei pa 100$ (n » 50), AE is approximately equal to 
10~ 3 <J) whereas it reaches 5 x 10~ 3 <I> for Ei w 200$. Nev- 
ertheless, we have verified that AE is smaller than the 
energy level spacing indicating that PT remains valid, 
even for quite large energies. As shown in the inset, 
the curve AE vs Ei behaves almost periodically, with 
a period approximately equal to ft, indicating that the 
PT accuracy depends on the nature of the unperturbed 
states. PT is exact for the unperturbed states that in- 
volve the exciton state insensitive to the phonons. By 
contrast, unperturbed states located in the neighborhood 
of quasi-degenerate states are the less well-corrected. 

As shown in Fig. 4, AE slightly increases with x 
so that PT remains valid in the intermediate coupling 
regime. It reproduces the influence of the exciton-phonon 
coupling on the energy levels up to x = 20 pN (77 = 1.41 
cm -1 ). In particular, PT accounts for the energy level 
crossing process that affects the dressed states. Owing to 
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the coupling V, unperturbed states hybridize giving rise 
to energy shift and splitting characteristic of anti-crossing 
phenomena. Some energy levels repel each other whereas 
other energy levels get closer. These levels describing 
exact uncoupled eigenstates, energy level crossing takes 
place. 

The e dependence of the exciton energy correction Suj^ 
is shown in Fig. 5. Owing to the dressing by virtual 
phonons, the exciton energy is redshifted and Suj^ de- 
creases linearly with Eg. However, this behavior de- 
pends on the exciton state. For the state \tpo) insensitive 
to the phonons, 5u = 0. By contrast, the stationary 
waves fj, ^ o, ± experience an energy correction approx- 
imately equal to — 0.2Eb- This shift reduces to — O.IEb 
for the states |^o) and whose energies define the 

band edges, and for the quasi-degenerate states \ip±). 
Moreover, Fig. 5 reveals that the energy correction of 
the stationary waves is almost e independent. This is 
no longer the case for the quasi-degenerate states be- 
cause 5u>+ (resp. decreases (resp. increases) lin- 
early with e (not well distinguishable in Fig. 5). Note 
that Su} + = (5w_ = (HOSE'S when e = 0. 

Fig. 6 shows the e dependence of the difference be- 
tween the energy w M of a dressed state \xn) and the cor- 
rected energy oj m + Sw^ of a bare state For the sta- 
tionary waves /i^o,±,a correspondence occurs between 
both energies suggesting that the dressing mainly induces 
an energy renormalization without significantly modify- 
ing the state. In fact, we have verified that \xp,) ~ |Vv) 
for /i 7^ ±. Note that the correspondence is exact for 
the state insensitive to the phonons so that one defines 
a dressed state \ Xo ) = \ipo) { v = L/2 = o). By contrast, 
the correspondence disappears for the quasi-degenerate 
states and the smaller e is, the larger is the energy differ- 
ence. A strong hybridization occurs between \ip±) that 
results in the formation of two quasi-degenerate dressed 
states. Denoted \x±) { v = L/2 =F 1 = ±), these dressed 
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FIG. 6: (Color online) d) M — oj m — Suj^ vs e for X — 10 pN and 
L — 10. 
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FIG. 7: (Color online) e dependence of Cj± (full lines) and ui± 
(dashed lines) for x = 10 pN and L = 10. 



states mainly correspond to superimpositions involving 
\ip + ) and \ip_). 

The e dependence of the energy u>± of the quasi- 
degenerate dressed state |x±) is shown in Fig. 7. For 
quite large e values, Cj+ (resp. u)_) is slightly red-shifted 
when compared with w + (resp. In that case, 

Cj± follows the e dependence of the corrected bare en- 
ergy lo± + 5ui± and \x+) (resp. |%— )) basically cor- 
responds to (resp. \ip-))- As e decreases down 
to zero, decreases and it reaches ujq. By contrast, 
w_ increases and it converges to a value located below 
the band center. The energy difference characterizes an 
anti-crossing phenomena so that the energy levels of the 
dressed states repel each other. This is the signature 
of the hybridization between \ip+) and \ip~) that results 
in the formation of two quasi-degenerate dressed states 

| x ±)«(hM±hM)/V2. 

The e dependence of the phonon frequency correction 
(5f2„ is shown in Fig. 8. The phonon frequency is ei- 
ther red-shifted or blue-shifted depending on the nature 
of the exciton. A redshift is induced when the exciton 
occupies a dressed state |%„) with v = 0, L/2 — 1. The 



7 



o -1 

o 





v=10 




v = 9 




v = S 




v=7 




v = 6M 

v=5(o) 




v=4(+) 




v = 3 




v = 2 




v=l 




v = 

I I I I 



0.00 0.01 0.02 0.03 0.04 0.05 



E 

FIG. 8: (Color online) SQ. V vs e for x = 10 pN and L = 10. 



closer to the band edge the exciton energy is located, the 
larger is the shift. By contrast, an exciton that occupies 
a state \xu) with v — L/2 + 1,...,L yields a blueshift 
of the phonon frequency. As previously, excitonic states 
near the band edge favor the largest phonon frequency 
shift. Note that in the state \xo), the exciton produces 
no phonon frequency shift. As shown in Fig. 8, when 
the exciton occupies a state isomorphic to a stationary 
wave, the phonon frequency shift is e independent. By 
contrast, a strong e dependence occurs when the exciton 
is in a quasi-degenerate dressed state |x±). In that case, 
<50_(_ < increases when e decreases whereas <5f2_ > 
decreases when e decreases. For nonvanishing e values, 
it turns out that <5f2_ w — 5CI+. However, when e tends 
to zero, Sfl+ converges to zero whereas <5f2_ reaches a 
quite small nonvanishing value. In other words, when 
e — > 0, the frequency shift is negligible when the phonon 
is accompanied by an exciton in either \xo), \x+) or Ix-)- 

The key ingredients entering PT being characterized, 
let us now study the effective exciton propagator |Gi,o(i)| 
whose time evolution is displayed in Fig. 9. The fig- 
ure shows that PT perfectly agrees with exact calcu- 
lations over a long-time scale. Initially equal to zero, 
|G\trj(t)| first increases with time. Then, at time Tm, 
it reaches a maximum value Gm quite close to unity. 
Finally, it develops damped oscillations that fluctuate 
around 1/2. These oscillations support a high-frequency 
small-amplitude modulation whose behavior depends on 
e. When e = 0.020, |Glo(OI exhibits a double maximum. 
At t = 248.47$~ 1 , |Glo(*)| first reaches a local maximum 
whose value is equal to 0.83. Then, the absolute maxi- 
mum G M = 0.88 occurs at T M = 495.48$ _1 . Such a 
temporal structure is very sensitive to e. When e slightly 
varies around 0.020, the absolute maximum jumps be- 
tween the two values of the double maximum so that 
Tm exhibits a discontinuous character with respect to e. 
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FIG. 9: (Color online) Time evolution of |G i0 (t)| for x = 10 
pN, L = 10 and T = 300 K. Exact calculations (dashed lines) 
and PT (full lines). 
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FIG. 10: (Color online) Max[\G L o(t)\] vs e for x = 10 pN and 
L — 10. Exact calculations (dashed lines) and PT (full lines). 



By contrast, when e = 0.013, |Glo(*)| exhibits a single 
absolute maximum whose behavior remains quite stable 
when e varies. It occurs at Tm — 699.85$ _1 and its value 
Gm = 0.97 is very close to unity. 

The e dependence of the maximum value Gm = 
Max[|Gio(i)|] is illustrated in Fig. 10 for different tem- 
peratures. The figure shows that PT provides results in 
a quite good agreement with exact calculations. A small 
discrepancy occurs at high temperature and for small e 
values, only. When e decreases from 0.05, the curve Gm 
vs e exhibits a series of minima and maxima whose value 
depends on the temperature. A local minimum corre- 
sponds to a singularity. The curve exhibits a kind of 
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FIG. 11: (Color online) Max[\G L o{t)\] vs T for x = 10 pN 
and L — 10. The calculations have been carried out using 
PT. 



cusp in the neighborhood of which the time Tm (e) is dis- 
continuous. By contrast, a local maximum corresponds 
to a well defined point where the first derivative of Gm (e) 
vanishes. Close to a local maximum, Tm(c) remains con- 
tinuous. Whatever the temperature, it turns out that 
the absolute minimum of Gm(c) occurs for e ~ 0.021. 
Its value slightly decreases with the temperature and it 
varies from 0.86 for T = 100 K to 0.84 for T = 300 K. By 
contrast, the absolute maximum takes place for e rs 0.013 
V T. Equal to 0.99 for T = 100 K, it decreases to 0.97 for 
T = 300 K. Note that if one defines P = 100 X (1 - G M ) 
as the percentage of the lost information, one obtains 
P sa 3 % at high temperature. When e tends to zero, the 
maxima and the minima are becoming more frequent but 
less pronounced. Gm tends to a temperature dependent 
value approximately equal to 0.99 and 0.95 for T = 100 
and 300 K, respectively. We have verified that the curve 
Gm vs e depends on both the lattice size and the exciton- 
phonon coupling strength. These parameters modify the 
position and the value of the minima and of the maxima. 
For instance, for e = 0.01 and T — 300 K, one obtains 
Gm = 0.93 for L = 10 and x = 10 pN, G M = 0.89 for 
L = 20 and X = 10 pN and G M = 0.73 for L = 10 and 
X = 20 pN. Increasing either L or \ reduces the QST 
fidelity. 

Finally, the temperature dependence of Gm is illus- 
trated in Fig. 11 for different e values. The calcula- 
tions have been carried out using PT that is particularly 
suitable for the considered e values. When e = 0.021, 
the curve Gm(c) lies in the neighborhood of its absolute 
minimum. Consequently, over the temperature range dis- 
played in the figure, Gm takes the smallest values. Al- 
most uniform, it is approximately equal to 0.86 ± 0.01. 
The corresponding lost information is thus about 14 %. 
For e = 0.01 and 0.02, a critical temperature T* dis- 
criminates between two regimes. When T < T* , Gm 
rapidly decreases with the temperature. By contrast, 
when T > T*, smooth variations take place. Note that 



T* w 40 K for e = 0.01 whereas T* « 130 K for e = 0.02. 
Finally, for e = 0.013, the curve Gjvf(e) lies in the neigh- 
borhood of its absolute maximum. Consequently, Gm 
slightly decreases with the temperature. It varies from 
(almost) unity for T = 10 K to 0.97 for T = 300 K. It 
scales as G M ml— (T/T ) 2 with T Q w 1510 K so that 
the lost information remains smaller than 3 %. 



IV. DISCUSSION 

According to the numerical results, PT is a powerful 
tool for describing the exciton-phonon dynamics. Within 
the nonadiabatic weak-coupling limit, it accurately ac- 
counts for the spectral properties of the system over a 
broad energy scale. Moreover, PT is particularly suitable 
for characterizing the time evolution of the effective exci- 
ton propagator over a long-time scale, even at high tem- 
perature. In that context, both exact and approximate 
calculations have revealed the potential powerfulness of 
the proposed communication protocol. The key point 
is that its efficiency strongly depends on the coupling 
e between the QCs and the CC. When e is judiciously 
chosen, it turns out that the maximum value of |Glo(^)| 
is quite close to unity. A high-fidelity QST occurs over 
a broad temperature range, the lost information during 
the transfer remaining smaller than 3 %. By contrast, 
specific e values induce a hole in the transmission curves. 
The maximum value of the effective exciton propagator 
deviates from unity resulting in the impoverishment of 
the transfered information. The lost information drasti- 
cally increases and it can represent almost 15 % of the 
initially implemented infromation at high temperature. 
Of course, these results depend on both the lattice size 
and the exciton-phonon coupling strength whose the in- 
crease reduces the fidelity of the transfer. 

To interpret these results, PT can be used for deriving 
a simplified expression of Gi,o(t). Indeed, our numeri- 
cal studies have revealed that the so-called diagonal ap- 
proximation works quite-well. Consequently, the exact 
exciton-phonon eigenstates \^>i) basically reduce to the 
effective Hamiltonian eigenstates \x v ) <S> \ri). The trans- 
formation U in Eq. (flU)) behaves as the unit operator and 
the exciton-phonon entanglement mainly results from the 
modification of the exciton states {tp^) — > \Xu)> the renor- 
malization of the exciton energies uj^ — > w„ and the cor- 
rection of the phonon frequency £1 — >• ft + £fi„. GLo(i) is 
thus rewritten as 



G L0 (t) m^Fuity-^mxuKx.lO), 



(14) 



where F u (t) — Zg\t)/ZB is the decoherence factor as 



F v {t) 



1 



,-pn 



(15) 



In Eq. (fT4"|) . Glo(£) is the sum of the probability ampli- 
tudes associated to the different paths that the exciton 
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can follow to tunnel between the QCs. A given path de- 
fines a transition through the dressed state \xv)- The 
corresponding amplitude involves the weight of the lo- 
calized states in the dressed state, a phase factor that 
accounts for the free evolution of the dressed exciton and 
the decoherence factor. This later contribution originates 
in the exciton- phonon interaction. Indeed, when Eq.® 
is developed in the phonon number state basis, GLo(t) 
can be interpreted as follows. The system being prepared 
in the factorized state |0) ® \n), Glo(*) is the probabil- 
ity amplitude to observe the system in a factorized state 
\L) ® exp(— inClt)\n). It thus describes an excitonic tran- 
sition during which the phonons evolve freely. When the 
exciton occupies the state |x„), the phonon frequency is 
modified. The probability amplitude that the phonons 
evolve freely reduces to a phase factor exp(— in8VL v t). 
Of course, such a phase factor does not affect the exci- 
tonic coherence when the phonons are initially in a pure 
state. However, at finite temperature, the phonons are 
described by a statistical mixture of number states. A 
thermal average is required and it yields a sum over dif- 
ferent phase factors which interfere with the other. The 
interferences lead to a decay of the excitonic coherence 
encoded in the decoherence factor. 

At this step, Eq.([I3)l can be simplified because only 
the dressed states that involve the localized states |0) and 
\L) contribute significantly to the exciton propagator. In- 
deed, the bare states that correspond to the station- 
ary waves of the CC (/x 7^ o, ±) are not degenerated. The 
dressing by virtual phonons mainly induces energy renor- 
malization without significantly modifying their nature. 
Therefore, the system exhibits N — 1 dressed states iso- 
morphic to these stationary waves, i.e. \Xfi) ~ an d 
d)p ps ujfj, + <5w p for jU 7^ o, ±. Because these dressed states 
are almost independent on the localized states, their con- 
tribution to Gio(i) can be disregarded. Note that |^) is 
coupled with z^ neighboring states through phonon ex- 
changes. In the nonadiabatic limit (\ujp, — lu^±i\ <§C O), 
Eq.(fTUl) yields Suj^ ps —z^Eb/L, as observed in Fig. 5. 
Because z% = Zl = 1, the shift experienced by the band 
edge states fi = 1 and \i = L is two times smaller than 
the shift experienced by the remaining stationary waves 
for which z^ = 2. In a market contrast, the bare state 
\if) ) does not interact with the phonons. The correspond- 
ing dressed state \x ) = IV'o) only involves the localized 
states |0) and \L) (see Eq.©) and the decoherence fac- 
tor F {t) reduces to unity. Consequently, \x ) defines an 
ideal path for the excitonic transition between the QCs. 

The fundamental point concerns the quasi-degenerate 
states \ip±) that are profoundly perturbed by the dress- 
ing mechanism. As for the stationary waves, the exciton- 
phonon interaction induces a redshift 6ui± of each energy 
ui±. This interaction originates in the dependence of \ip±) 
with respect to the robust stationary wave \ <Pl/2) ■ There- 
fore, each quasi-degenerate state is coupled with two sta- 
tionary waves through phonon exchanges. In that case, 
it is easy to extract Sui± from Eq. fTUl) . In doing so, the 
shifts can be expanded as a Taylor series with respect to 



e as 5v± = — YlrLo(~ ^) r E r e r . The positive coefficients 
E r are defined as 



Er 



E r 



E r 



2 \(n + Auj) r + 1 (Q-Aui) 



r+l 



(16) 



J L/2 



W L/2±ll' In the 



where E = 2$/\/L and Alu = 
nonadiabatic limit, Eq ps —Eb/L so that 8ui+ (resp. 
<5cj_) decreases (resp. increases) linearly with e form 
—Eb/L, in a quite good agreement with the results dis- 
played in Fig. 5. Note that 6lj± is similar to the shift 
experienced by the band edge states because the weight 
of \v>l/2) m IV'i) is equal to \/\[2. 

In addition to energy renormalization, the exciton- 

phonon interaction yields an effective coupling = 

-(Soj + + Suj-)/2 between \ip+) and |^_) (Eq.JTD])). Quite 
similar in magnitude to the energy shifts, this coupling 
can no longer be neglected resulting in a strong hybridiza- 
tion between the quasi-degenerate states. These states 
behave as a two-level system independent on the remain- 
ing states whose eigenstatcs define two quasi-degenerate 
dressed states as 



| X +) w +cos(0)|^ + )+sin(W_> 
|X_> w -sin(fl)|V+>+cos(fl)|^_). 

The corresponding eigenenergies are expressed as 

Sui + + <5w_ 1 r— 
uj ± = uj H ± 2^ A ' 



(17) 



(18) 



where 



A = (w+ + Suj + - u- - 5uj-) 2 + Av 2 + _ 



cos(26») 



UJ + + 5iJ+ — LJ- — SiO- 



(19) 



The hybridization is encoded in the 9 parameter whose 
value depends on both the exciton-phonon coupling 
strength, measured by 5u>±, and the energy difference 
S = u> + — ui- = 2Ee. Two asymptotic situations oc- 
cur. For quite large e values, \Suj±\ < i so that the 
influence of the quasi-degeneracy is negligible. The hy- 
bridization is weak and the dressed states reduce to the 
bare states, i.e. \x±) ~ \4 > ±) arL d £>± ~ u>± + 5w±. 
The bare energies experience a similar redshift approx- 
imately equal to —Eb/L, as displayed in Fig. 7 for 
e ps 0.05. By contrast, for small e values, \Su>± \ ^> 5 and 
the hybridization is strong. One obtains 9 ps 7r/4 and 
\x±) ps (1^+) ± |V'-))/v / 2- The energy _w+ reaches the 
band center and it scales as w+ ps ujq + (E — Ei) 2 e 2 /2Eq. 
By contrast, w_ ps L}Q-2E Q -2E 2 e 2 -(E-Ei) 2 e 2 /2E lies 
close to ujq — 2Eb/L, in a good agreement with the results 
displayed in Fig. 7. Note that in that case |x+) reduces 
to a superimposition of the localized states whereas |x_) 
tends to the robust stationary wave. 

By combining Eqs.© and (TT71) . it is easy to show that 
the quasi-degenerate dressed states depend on the local- 
ized states. They thus define two relevant paths for the 



10 



excitonic transition between the QCs. However, when 
the exciton occupies |x±), it induces a shift 5il± of the 
phonon frequency. Quantum decoherence is no longer 
negligible and the decoherence factors is expressed as 



-in8Q+t 



F±(t) 



1 + 4 An 2 sm 2 (Sn±t/2) 



(20) 



where n — [exp(/3f2) — 1] _1 and An 2 = h(fi + 1). Note 
that Eq. ([20"l) has been obtained from Eq. ([T5|) by simpli- 
fying the time dependence of the argument of the deco- 
herence factor. By inserting Eq. (fT7l) into Eq. flTUl) . the 
phonon frequency shift is written as 



Sn± « =p2 cos(26>) E 2r +ie 



2r+l 



(21) 



r=0 



In a quite good agreement with the results displayed in 
Fig. 8, Eq. (l2~Tj) reveals that when the exciton occupies 
the state |x+) (resp. |x-))j it induces a redshift (resp. 
blueshift) of the phonon frequency. Moreover, as ob- 
served in Fig. 8 when e is not too small, one obtains 
<5f2_ = — Sn + . However, when e tends to zero, Eq. (|2"Tj) 
yields 8Q,± s» ^f2E 1 (E - E 1 )e 2 /E Q . Such a behavior 
slightly differs from our numerical results. As shown in 
Fig. 8, <5i7+ vanishes when e tends to zero whereas <557_ 
reaches a rather small value. This discrepancy originates 
in the fact that the exact dressed states |x±) slightly de- 
pend on the stationary waves \ipn), with /i ^ o, ±. When 
e tends to zero, such a dependence is more pronounced for 
|x_). These additional components, that are not taken 
into account in Eq. (|17|) . yield a nonvanishing <5f2_ value. 

In that context, it turns out that |xo), |x+) an d |X-) 
define the main paths followed by the exciton to tunnel 
between the QCs. Consequently, Gi,o{t) is rewritten as 



Gzo(t) ~ -Ajvcxp(-iw i) 



(22) 



-\F + (t)\exp(-iW+t)[l + sm(2e)] 
-|F_(t)|exp(+iW_t)[l-sin(20)] 



where W± = ±(w± — ujq + fi5Q±) defines a positive fre- 
quency relative to the band center lu q . According to 
Eq.([2"2"j). (j\Lo(t) is the sum of the three transition am- 
plitudes, t = 1/2 and r± = -\F±{t)\exp(^iW±t)[\ ± 
sin(20)]/4, connected to the three relevant paths. Be- 
cause |xo) is insensitive to the phonons, the probability 
amplitude r Q is the weight 1/2 of the localized states. By 
contrast, the probability amplitude r± that the exciton 
tunnels through |x±) involves a phase factor whose time 
evolution is governed by T± — tt/W±. This phase fac- 
tor is weighted by both the decoherence factor and the 
weight of the localized states. Consequently, depending 
on the value of the model parameters, W±, and F±(t) 
will take particular values so that each probability am- 
plitude will develop a specific time evolution. Quantum 
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FIG. 12: (a) G M (e) vs e for X = 10 pN, L = 10 and T = 300 
K using PT. (b) a = W s /W f vs e. (c) T M /T f vs e. 



interferences between the different amplitudes will occur 
resulting in a characteristic time evolution of |G\i,o(i)|- 
As time evolves, |G/,o(£)| will reach a maximum whose 
value Gm will depend on the model parameters through 
the interference pattern. 

To understand this phenomenon intuitively, let us fo- 
cus our attention on the influence of the coupling e be- 
tween the QCs and the CC. Provided that e is not too 
small, the numerical results suggest that <55!_ m — 5Q+ > 
0. The modulus of the decoherence factor is thus al- 
most independent on the nature of the exciton so that 
\F v (t)\ w F(t) Vis = ±. Moreover, it turns out that the 
influence of 6 remains negligible, i.e. sin(20) f=s 0. There- 
fore, after straightforward algebraic manipulations, one 
obtains 



1, 



F 2 (t) cos 2 (W a t) (23) 
2F{t) cos{W s t) cos(W>t)], 



where W s = (VF_ - W + )/2 and W f = (TF_ + W + )/2. In 
the right hand side of Eq. (|2"31 , the first term is the prob- 
ability |t | 2 that the exciton realizes a transition through 
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|Xo). The second term reduces to |r + + r_| 2 and it refers 
to the probability that the exciton tunnels through either 
\x+) or |x-)- The last term describes the quantum in- 
terferences that arise when the exciton follows the paths 
involving either |xo) or \x+), and |xo) or |x-}- It mixes 
W + and W- and it behaves as an amplitude-modulated 
signal f(t) = F(t) cos(W s t) cos(Wft). It thus exhibits a 
slowly varying envelope with frequency W s that supports 
a fast modulation with frequency Wf, both components 
being damped by the decoherence factor F(t). Such a be- 
havior allows us us to introduce the fundamental times 
T f = ir/Wf and T s = n/W s . 

The time evolution of |Gio(i)| is mainly encoded in the 
quantum interference term f(t). Consequently, |G io(*) I 
reaches its maximum value Gm(c) for a time Tm(e) that 
minimizes f(t). As illustrated in Fig. 12(c), such a situ- 
ation occurs when Tm is a multiple of Tf, i.e. Tm = nTf 
with n — 1,2,.... However, the integer n and the mag- 
nitude of Gm(c) depend on e through the value of the 
ratio a = W s /Wf. This ratio extends from zero for large 
e values to unity for small e values, as displayed in Fig. 
12(b). Therefore, as e varies, different regimes take place. 

For specific e values, a is such that a situation arises 
in which the quantum interference term f(t) exhibits two 
equivalent minima at two distinct times. In that case, 
the curve Gm vs e exhibits a local minimum and the 
time Tm shows a discontinuity between two multiples 
of Tf. We have observed that in that case there are 
two equivalent configurations for which the probability 
amplitudes r + and r_ are in phase with each other but 
are not exactly in phase with the amplitude t . At high 
temperature, such a situation occurs for the e values that 
provide a = (2q- l)/(2q + 1) with q = l,2, .... The time 
Tm exhibits discontinuities between qTf and (q + 1)2/ 
and one obtains T + = 2gT_ and T f = 2T + /(2q + 1). 
Note that as shown in Fig. 12, the largest hole in the 
transmission curve occurs for q — 1, i.e. for a = 1/3, so 
that Tm varies between Tf and 2Tf. 

In a marked contrast, for particular e values, a is such 
that the quantum interference term f(t) exhibits an iso- 
lated minimum close to —1. Such a situation takes place 
at a time Tm(c) so that the probability amplitudes r+ 
and r_ are simultaneously in phase with the probability 
amplitude r Q . Constructive interferences occur between 
the different paths followed by the exciton to tunnel be- 
tween the QCs. As a result, the curve Gm vs e exhibits a 
local maximum (see Figs. 10 and 12(a)). Such a situation 
appears for the e values that yield a = {q — p)/{q+p+ 1), 
where p and q are two positive integers such that q > p. 
One thus obtains T M = (2p + 1)2+ = {2q + 1)T_, i.e. 
T M = {p + q + 1)2/. In that case, f(T M ) = -F{T M ) so 
that the maximum value of the exciton effective propa- 
gator is approximately equal to Gm = [1 + F{T M )]/2. 

As shown in Figs. 10 and 12(a), the largest local max- 
imum occurs for a specific value e*. When e = e* , p = 
and q = 1 so that a = 1/2 and T M = T + = 3T_. Solv- 
ing the equation T + = 3T_ by exanding the parameters 
with respect to e provides an analytical expression for e* 



written as 



V2E 



^{E - E x ) 2 - AnEi(E - E x ) - 2E a E 2 



(24) 



Over the temperature range T = 100 - 300 K and for 
X = 10 pN and L = 10, Eq.^ yields e* 0.011 
in a quite good agreement with the numerical estimate 
e* sa 0.013 (Fig. 10). Therefore, when e = e*, T+(e*) 
defines the shortest time for which constructive inter- 
ferences occur. The decay provided by the decoherence 
factor is thus minimized and the QST is optimized. In 
the nonadiabatic weak-coupling limit one thus obtains 



G M (e*) 



1 



tf 2 An 2 



(25) 



where 5n = uq — n and uq = \uj + (e*) — Wo|/|5f2+(e*)|. 
Eq. (f25"|) clearly shows that the impoverishment of the 
transfered information results from the thermal fluctu- 
ations An of the phonon number. Fortunately, in the 
weak coupling limit, no 3> n so that An is always smaller 
than Sn. Because An is proportional to hsT/Sl, the op- 
timized value of the effective exciton propagator scales 
as G M (e*) «1- (T/T ) 2 with T = 2n f2/7r. As ob- 
served in Fig. 11, Gj\/(e*) only slightly deviates from 
unity provided that T remains smaller than the critical 
temperature To. With the parameters used in the simu- 
lation (x = 10 pN and L = 10), T m 1600 K indicating 
that the lost information during the transfer is negligible, 
even at room temperature. 

Consequently, when e is judiciously chosen, construc- 
tive interferences take place between the different paths 
followed by the exciton to tunnel between the QCs. The 
influence of the quantum decoherence is minimized and 
an ideal QST occurs over a broad temperature range. 
The fidelity of the transfer remains quite close to unity, 
even at high temperature, indicating the powerfulness of 
the proposed communication protocol. 



V. CONCLUSION 

In the present paper, a new communication protocol 
has been proposed in which QST is achieved by a high- 
frequency vibrational exciton. The main idea was to use 
two distant molecular groups grafted on the sides of a 
molecular lattice. These groups behave as two QCs on 
which the information is encoded and received. By con- 
trast, the lattice defines the CC along which the exciton 
propagates and interacts with a phonon bath. To mini- 
mize the impact of the quantum decoherence, the struc- 
ture was designed so that a vibrational resonance occurs 
between the QCs and the robust stationary wave of the 
lattice whose energy is exactly located at the band center. 

To highlight the relevance of PT, special attention has 
been paid for describing a simple model in which an 
exciton is dressed by a single phonon mode, only. In 
that case, the Hamiltonian was solved exactly so that 
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the PT accuracy has been checked. Within the nonadi- 
abatic weak-coupling limit, it has been shown that PT 
is a powerful tool for characterizing the exciton-phonon 
dynamics. Provided that the lattice size is not too large, 
it yields a very good estimate of the spectral properties 
over a broad energy scale. Moreover, it has been observed 
that PT is particularly suitable for describing the time 
evolution of the exciton RDM, even at high temperature. 

In that context, it has been shown that the system 
supports three quasi-degenerate exciton states that in- 
volve the states localized on the QCs. When the exciton 
occupies one of these states, it does not significantly mod- 
ify the phonon bath and keeps its coherent nature over 
a long-time scale. These states define the relevant paths 
followed by the exciton to tunnel between the QCs. Con- 
sequently, when the coupling between the QCs and the 
CC is judiciously chosen, constructive interferences take 
place between these paths. The quantum decoherence 
is minimized and an almost ideal QST occurs. The fi- 
delity of the transfer remains quite close to unity over a 
broad temperature range, indicating the powerfulness of 
the proposed communication protocol. 

Finally, because the present approach has clearly re- 
vealed the relevance of PT, it will be generalized for de- 
scribing the influence of all the phonon modes. The main 
idea is to check the efficiency of the proposed protocol 
with a more realistic system whose dynamics cannot be 
solved exactly. 

Appendix A: Quasi-degenerate second order 
perturbation theory 

Quasi-degenerate PT involves a unitary transforma- 
tion U — exp(S) that provides a block-diagonal trans- 
formed Hamiltonian H = UHU< in the unperturbed ba- 
sis „). More precisely, the desired Hamiltonian must 
be diagonal with respect to the phonon number states, 
only. To proceed, any operator O acting in £ can be split 
as O = Od + OnDi where Od is the diagonal part with 
respect to the phonon number states whereas On d is the 
remaining non-diagonal part. Note that such a partition 
is equivalent to that provided by a projector formalism 59 . 
In the unperturbed basis, these operators are defined as 

(K,n\OnK>, nt ) = (K,n\°K',n)8nn> (Al) 

In that context, because Vd = 0, one seeks the anti- 
hermitian generator S = Snd as a non-diagonal operator 
with respect to the phonon number states. It is expanded 
as a Taylor series as S = Si + S2 + ••■ where S q is the 
qth order correction in the coupling V. Consequently, H 
becomes 

H = H + V + [Si,H ] 

+ [S 1> V} + [S 2 ,H } + ~[S 1 ,[S 1> H }} + ... (A2) 



From Eq. (|A2[) . S is derived order by order to obtain a 
block-diagonal form for H at the desired order. Up to 
second order, the solution is given by the equations 

[Hq, Si] = Vnd 

[#0,S 2 ] = - [S\,V]nd 

H = H + ±[S 1 ,V] D . (A3) 

Because V — M(a' 1 + a) is a linear combination of 
creation and annihilation phonon operators, Si is of the 
form Si = Za 1 ' — Z^a. The unknown operator Z acts in 
£a, only. No restriction affects this operator because 
S{ = —Si- Therefore, inserting this expression into 
Eq.{S9 yields Z^ = M^/fa/i - + n )- 

The knowledge of Si allows us to compute the commu- 
tator [Si,V"] that is required to derive both H and S2. 
This commutator is defined as 

- [Si, V] = A + Ba)a) + B^aa + (B + B*)a!a, (A4) 

where B = [Z,M]/2 and A = -{Z^M + MZ)/2. From 
the diagonal part of Eq. (|A4j) . H reduces to 

H = H A + A + (Q + B + B r )a <! a. (A5) 

We thus recover Eq.© with SH A = A and <5£1 = B + 
B' whose representation in the unperturbed basis yields 
Eq. (fTU|) . From the non diagonal part of Eq. (|A4[) . one 
seeks S2 of the form S2 = Ea^a^ — E^aa. The unknown 
operator E acts in £a, only. Inserting this expression 
into Eq.(|A3} yields E^, = - uy + 20). 



Appendix B: Effective exciton propagator 

Expanding U as a Taylor series with respect to V, 
Eq. (fT"3]) yields the second order expression of the effec- 
tive exciton propagator as 

G L0 (t) = £^W e -^x 

+ (L\Z\x„)(x»\Z i \0)n^(ty( n+m » )t 

+ (L\Z^\ Xv )( X u\Z\0)(n^(t) + i) e -*Cn+«>,)t 

- (L\Z&\x»){Xv\0)nW(t)/2 

- (L\Z^Z\ X u){Xu\0){n^(t) + l)/2 

- {L\xu)<X„\ZZl\0)nM{t)/2 

- (L\xu){Xu\Z^Z\0){n^(t) + l)/2 }, (Bl) 

where n^\t) = [cxp(/30 + i8tl„t) - l] -1 . 
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